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ABSTRACT 


A laboratory simulation experiment of the interaction of the solar wind and the 
earth’s magnetic field is being conducted at Lewis Research Center of NASA. The 
plasma flow is produced by an MPD arc and a magnetic dipole is used to simulate the 
earth’s magnetic field. The initial investigations of the stagnation region of the flow 
show a rapid rise in the plasma potential. The magnetic field in this region is also al- 
tered from the dipole type behavior it exhibits when no flow exists. This region of 
change in plasma potential and magnetic field has a thickness characterized by the elec- 
tron cyclotron radius. It appears to be a merging together of a collisionless decelera- 
tion and stagnation of the flow. However, whether or not a collisionless shock wave 
exists has not been conclusively demonstrated. 
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SUMMARY 

A laboratory simulation experiment of the interaction of the solar wind and the 
earth’s magnetic field is being conducted at Lewis Research Center of NASA. The plas- 
ma flow is produced by an MPD arc and a magnetic dipole is used to simulate the earth’s 
magnetic field. The initial investigations of the stagnation region of the flow show a 
rapid rise in the plasma potential. The magnetic field in this region is also altered from 
the dipole type behavior it exhibits when no flow exists. This region of change in plasma 
potential and magnetic field has a thickness characterized by the electron cyclotron ra- 
dius. It appears to be a merging together of a collisionless deceleration and stagnation 
of the flow. However, whether or not a collisionless shock wave exists has not been con- 
clusively demonstrated. 


INTRODUCTION 

For the past few years NASA LeRC has been conducting preliminary experiments on 
the laboratory simulation of the interaction of the solar wind with the earth's magnetic 
field (ref. 1). Recently a detailed plasma diagnostics study of these experiments has 
been initiated. This report will report on the results of this diagnostics program. 

The electric thruster vacuum facilities available at Lewis are ideal for such an ex- 
periment. A large 15-foot diameter by 6 5- foot long vacuum tank houses the experiment. 
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The pumping speed of this tank is 500 cubic meters per second in the 10 - to 10 -torr 

pressure range. 

To simulate the solar wind a magnetoplasmadynamic (MPD) arc was used. Hydro- 
gen is the propellant gas. This device is being studied as a possible thruster for deep 
space missions (ref, 2). A magnetic dipole is used to simulate the earth’s magnetic 
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field. The axis of the dipole is alined perpendicular to the plasma flow. 

Initial investigations have been confined to the stagnation region of the flow where a 
collisionless shock wave, similar to the earth’s bow shock wave, should exist. The 
structure of this region was determined by measuring the plasma potential and magnetic 
field profiles. Currently, the electron number density and temperature profiles are also 
being investigated. 


DIAGNOSTICS 

Figure 1 shows a side view of the experimental apparatus. The magnetic dipole is 



Figure l. - Diagnostic equipment for solar wind simulation experiment. 


enclosed in an insulating spherical shell. Shown are the three probes that are being used 
to determine the plasma properties. They are mounted on a traversing rake that is mov- 
able in the axial, transverse, and vertical directions. The emissive probe is used to 
measure the plasma potential. The Hall magnetic probe is water cooled and is alined to 
measure the magnetic field in the z-direction. A 0. 001- inch-diameter tungsten wire 
Langmuir probe is being used to determine the electron density and temperature. The 
thrust and hydrogen flow rate were measured as described in reference 2. 
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EXPERIMENTAL CONDITIONS 


The values of the significant plasma properties that exist in the plasma flow are 
given in table I for both the experiment and the solar wind. To produce a collisionless 
plasma flow such as exists in the solar wind the characteristic collision mean free paths 
must be large compared to the dimensions of the experiment. Other conditions that must 
be satisfied to simulate the entire earth solar wind interaction are discussed in refer- 
ence 1. In this experiment the distance between the MPD arc and the dipole (fig. 1) is 
slightly over 1 meter and the diameter of the beam is about 1 meter . The values of the 
electron- ion and electron-neutral mean free paths (table I) are an order of magnitude 


TABLE I. - PROPERTIES OF PLASMA FLOW 


Plasma property 

Solar wind 
experiment 

Solar wind 

Flow velocity, m/sec 

5xl0 4 

5x10® 

Electron number density, n^, cm" 3 

2.6X10 10 

3 

Electron temperature, kTe 1 /q, V 

7 

7 

Upstream magnetic field, B^, T 

3.5X10 -4 

4x10“ 9 

Upstream sonic Mach number, M g 

1. 5 

15.0 

Ratio of kinetic pressure to magnetic 
pressure, j8 = ^kTe^ j (b 3 ^2^ 

0.8 

0.6 

Electron- Ion mean free path, m 

10 

10 5 

Electron-Neutral mean free path, m 

10 


Ion- Neutral mean free path, m 

1 


Debye length, m 

10"' 5 

10 


greater than the dimensions of the experiment. However, the ion-neutral mean free path 
is about the same size as the dimensions of the experiment. This is one of the reasons 
this investigation has been confined to the region close to the MPD arc where the hydrogen 
ions are less likely to collide with a background gas atom. Within this first 1/2 meter 
of the arc lies the region where the plasma flow is being decelerated and compressed by 
impingement on the dipole magnetic field. If the flow velocity is high enough a collision- 
less shock will occur in this region. The mean exhaust velocity of the MPD arc, u^, was 
determined from the measured thrust and mass flow. 

It is important to note that for this experiment the ratio of the kinetic pressure to 
the magnetic pressure, /3 = 0. 8, is close to that for the earth's bow shock, which is about 
0. 6. In the experiment of Pugh and Patrick (ref. 3) j3 has a much lower value of 0. 2. 
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RESULTS AND DISCUSSION 


Figure 2 shows the experiment in operation. The view is at an angle of approxi- 
mately 45° from slightly behind the MPD are. All the initial measurements have been 
made on the center line of the flow in the bright region just downstream of the MPD arc. 



Figure 2. - Solar wind simulation experiment. 


The plasma potential, V , relative to the tank-wall obtained with the emissive probe 
is shown in figure 3 as a function of distance from the dipole, x. In figure 4 the magnetic 
field perpendicular to the flow is presented as a function of x. The output of the emis- 
sive probe and the Hall magnetic probe were read on the y-axis of an x-y plotter and the 
position of the probe on the x-axis. As a result a continuous plot of V p and B z as a 
function of x is obtained. Two different profiles of V p are shown in figure 3. The only 
major difference between the two is that the region of rapid increase moves back about 
4 centimeters for the case where the dipole current is 495 amperes. This is as expected, 
since the stagnation magnetic field will be closer to the dipole for the lower current. 

Both potential profiles show a rise of about 12 volts in the initial discontinuity. This 
rising potential region may be a normal shock with the potential increase corresponding 
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Figure 3. - Plasma potential profile. 



Figure 4. - Magnetic field perpendicular to plasma flow. 
Dipole current 1000 amperes. 


to the velocity decrease associated with the shock. If this is so then 

i m o( u f - u !)= “ iV p 


a) 


where Uj is the upstream velocity, is the downstream velocity, m Q the ion mass, 
q the charge, and AV^ the change in potential. If AV^ = 12V, u^ = 5x1 0^ meters per 
second, then for hydrogen ions, equation (1) yields the following for Ug. 
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and 


u 2 = 4. 5xl0 3 m/sec 


The limiting velocity ratio for a normal MHD shock (ref. 4) is 4. Therefore, this poten- 
tial jump cannot represent a shock. However, it may represent the total velocity change 
that will occur across a shock and stagnation region together. The potential change re- 
quired to stagnate hydrogen ions moving with a velocity of 5x1 0^ meters per second is 
13. 1 volts. 

To assure the existence of a shock the upstream Mach number, Mp must be greater 
than one, where is defined as follows: 


M 


1 " 



where 


a 


1 



( 2 ) 


(3a) 


is the upstream sonic speed and 


V 


a 


B, 


V Mm o n l 


(3b) 


is the upstream Alfven speed. 

If a Mach number is calculated using the measured values of a^, Tp n^, and 
given in table I, then we find M^ ~ 0.9. However, if B^ = 2 gauss, which, according to 
figure 4, occurs slightly upstream of the discontinuity, then M^ ~ 1.2 From these re- 
sults it is not clear that M^ > 1 at the discontinuity. It may be that the flow is being 
compressed without the formation of a shock. 
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The thickness of the potential jump region is about 10 centimeters. The upstream 
electron cyclotron radius, 




( 4 ) 


while the ion cyclotron radius, 


u„ 


Vr 


U. 


qB. 


200 cm 


m i 


(5) 


The so-called Chapman-Ferraro-Rosenbluth collisionless sheath thickness (ref. 5), 
which is proportional to the ratio of the speed of light to the electron plasma frequency is 
the following. 



This quantity is the characteristic length for a beam plasma to be stagnated by a mag- 
netic field. Referring to equations (4) to (6), it appears that the electron cyclotron 
radius is the dimension that characterizes the thickness of this potential jump. This re- 
sult agrees with the model Tidman (ref. 6) presents for the earth’s bow shock. 

Pugh and Patrick (ref. 3) find that the characteristic dimension of their shock front 
is the ion cyclotron radius defined as follows. 


PP 



(7a) 


where, 



B„ 


V^ m o n l 
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is the Alfven speed. Equation (7a) therefore becomes 


PP 


c 

" °i/ V m e 


5 cfr ~ 87 cm 


(7b) 


where is the upstream ion plasma frequency. For this experiment r is much 

PP 

larger than the thickness of the discontinuity. 

Another interesting comparison with the Patrick and Pugh experiment (ref. 7) is 
that they find no potential jump through their shock front. Rather they find that the rela- 
tion, 


E + uXB = 0 


( 8 ) 


is satisfied, 
equation (8), 


So the streamlines of the flow represent equal potential lines. 



Also, from 


(9) 


The electric field, E^, was measured with the emissive probe and found to be E y ~ 0. 3 
volts per centimeter. Therefore for Bj ~ 3 gauss, 

u ~ 10 5 m/sec 

which is a factor of 2 larger than the measured value of u^. It therefore appears that 
equation (8) does not apply to this experiment. 

Figure 4 shows that when the thruster is not operating the magnetic field has the 

O 

usual 1/x dependence of a dipole field. With the thruster operating the magnetic field 
shows a change occurring in the same location as the rapid rise in the plasma potential. 
If we consider the region of change as lying between x ~ 69 centimeters and x ~ 59 cen- 
timeters then, 


B 2 8 . 0 _ 2 3 

B x 3. 5 

If the change in magnetic field is assumed to occur across a normal MHD shock, then 
from shock relations (ref. 4), this corresponds to ~ 2.2. This Mach number is 
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higher than expected. Thus the magnitude of the magnetic field change also supports the 
contention that this region of charge is not just a shock. 

Besides the fact that it is not known for certain that >1, the other significant 
point that arises in the experiment is that there does not appear to be a separation be- 
tween the shock front and the stagnation region. This point has already been discussed 
in connection with the plasma potential and magnetic field profiles. The shock, if there 
is one, is probably merged into the stagnation region. The rise in potential at x = 48 in 
the plasma potential profile (fig. 3) is probably behind the stagnation point since it occurs 
in the dark region of the flow (fig. 2). The magnetic field profile (fig. 4) shows no sec- 
ond ’’bump” that might indicate a stagnation point. Based on supersonic flow theory, the 
lack of separation between the shock and stagnation point is not surprising. For a sphere 
with M_ - 1.5, the detachment distance for the bow shock is about 1/2 the radius of the 
sphere (ref. 8). In this experiment the equivalent sphere radius is about 60 centimeters. 
Therefore, the detachment distance, based on supersonic flow theory, should be about 
30 centimeters. The region of rapid change has a thickness of 10 centimeters. There- 
fore, the condition that the detachment distance be much larger than the shock thickness 
in order to insure sepa ration between the shock and stagnation point is not satisfied. 

Some idea of the magnitude of the stagnation point magnetic field can be obtained by 
equating the stagnation point magnetic pressure to the free stream momentum flux plus 
pressure. 


B s 2 B 1 

_.m A . 1+h H 1+ _ 


Using the values given in table I, 


(10) 


B 0 ~ 7 gauss 
s 

Referring to figure 4 we see that the 7-gauss point occurs at the end of the discontinuity. 
Based on this result it appears that the stagnation region occurs within the region where 
the magnetic field deviates from the dipole field. 


CONCLUSION 

A region of rapid change in the plasma potential and magnetic field of a hydrogen 
plasma impinging on a magnetic dipole has been measured. The thickness of this region 
(10 cm) is characterized by the electron cyclotron radius. It appears that this region 
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j may be the merging together of a collisionless shock and the stagnation region of the 
j flow. Whether or not a shock exists, however, has not been conclusively demonstrated. 

Lewis Research Center, 

National Aeronautics and Space Adminisrration, 

Cleveland, Ohio, December 10, 1968, 

129-02-08-03-22. 
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